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Abstract. The large and saturating differential elliptic flow V2 (p±) observed in 
Au + Au reactions at RH1C so far could only be explained assuming an order of 
magnitude denser initial parton system than estimated from perturbative QCD. 
Hadronization via parton coalescence can resolve this "opacity puzzle" because it 
enhances hadron elliptic flow at large p± relative to that of partons at the same 
transverse momentum. An experimentally testable consequence of the coalescence 
scenario is that V2 (p± ) saturates at about 50% higher values for baryons than for 
mesons. In addition, if strange quarks have weaker flow than light quarks, hadron 
i>2 at high p± decreases with relative strangeness content. 



PACS numbers: 12.38.Mh; 24.85.+p; 25.75.Gz; 25.75.-q 

1. Introduction and summary 

Differential elliptic flow, V2(p±) = (cos(2 4>)) p± , the second Fourier moment of the 
azimuthal momentum distribution for a given p± , is one of the important experimental 
probes of collective dynamics in A + A reactions [Q. Measurements of elliptic flow at 
high transverse momentum provide important constraints about the opacity of the 
partonic medium in heavy- ion collisions [21 ED El- 

Recent data from RHIC for Au+Au reactions at ^/snn = 130 and 200 GcV show 
a remarkable saturation property of elliptic flow for ~ 2 GeV < p± <~ 6 GeV with 
V2 reaching up to 0.2 |51 |SJ |SJ [7] , a value that corresponds to a factor of two azimuthal 
angle asymmetry of particle production relative to the reaction plane. As illustrated in 
Fig-d the saturation pattern differs qualitatively from calculations based on a variety 
of different models: ideal (nondissipative) hydrodynamics"""] 111! 1101 I12| . hadronic 
cascades with string dynamics"""]""), inelastic parton energy loss [5], and classical Yang- 
Mills evolutional]. 

On the other hand, covariant parton transport theory [T51 1161 "T"*"l 1171 151 119 | is able 
to reproduce the observed flow saturation pattern. The main ingredient of this theory 
is a finite local mean free path A(s, x) = l/a(s)n(x), providing a natural interpolation 
between free streaming (A = 00) and ideal hydrodynamics (A = 0). Several studies 
show that initial parton densities and elastic 2^2 parton cross sections estimated 
from perturbative QCD, dN g /dr](b = 0) ~ 1000 ["""""j and Qgg^gg ~ 3 mb, generate too 
small collective effects at RHIC fTol 1171 15*1 IT5 j . Nevertheless, quantitative agreement 
with the «2(rl) data is possible, provided initial parton densities and/or cross sections 
are enhanced by an order of magnitude to adN g /di](b = Q) ~ 45 000 mb 3 . An 
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Figure 1. Illustrative comparison between V2(p±) data from RHIC and model 
calculations |2] [T21 IS1 ITT1 [1771 11211131 . 

enhancement of parton opacities seems necessary to explain the pion interferometry 
data from RHIC as well|T5|. 

To compare to the experiments, parton transport models have to incorporate 
the hadronization process. The studies mentioned above considered two simple 
schemes: 1) lparton — ► Itt hadronization, motivated by parton-hadron duality, and 2) 
independent fragmentation. An alternative is parton coalescence |21j. which became 
the focus of recent theoretical interest [22123 ESI S El as a possible explanation for the 
anomalous meson/baryon ratio and features of the elliptic flow data at RHIC. These 
studies indicate that hadron production at RHIC may be dominated by coalescence 
out to surprisingly large transverse momenta p± ~ 5 GeV. 

In this talk I show that hadronization via parton coalescence can resolve the 
opacity puzzle because it leads to an amplification of elliptic flow 4, at high p± (see 
Fig. EH)- With this enhancement, parton transport calculations 3 can account for the 
charged particle elliptic flow data from RHIC with moderate initial parton densities 
and cross sections dN g /drj ~ 1500 — 3000, <J gg ^gg = 3 mb. 

Parton coalescence predicts a rather simple relationship, Eq. J2J, between parton 
flow coefficients and those of hadrons. If all partons have similar elliptic flow, elliptic 
flow is stronger for baryons than for mesons, wf « 1.5u^, at high p± > 2 — 3 GeV, 
as shown in Fig. EH- If on the other hand strange quarks show weaker flow than light 
quarks, elliptic flow at high p± is ordered by relative strangeness content, such that 
v% > v£ ~ vf > v% > vtf > v$, > i>2 > v%, and > ~ 3«2/2 (see 

Fig- El 3 ) • These predictions are reinforced by preliminary calculations using the MPC 
parton transport model 26j (see Fig. EH) and can be readily tested in current and 
future heavy- ion collision experiments. 

With the help of (El , one can also determine the elliptic flow of the various parton 
species from V2 measurements for several hadrons, providing valuable insight on the 
evolution of the dense parton medium in heavy-ion collisions. 

2. Parton coalescence 



A convenient framework to start with is the local variant of the coalescence 
model 22, 4 , in which the coalescing partons have equal momenta. This is justified 
if hadron wave functions are narrow in momentum space and if constituent partons 
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have similar masses. See Ref. |23] for a study of the more general case. 

Assuming that all partons have the same momentum distribution, the hadron 
spectra at midrapidity are then given by 



|^(j?x) = C B (p ± ) 



dN q 



1 3 



2 



d 2 pj 

where Cm and C b are the probabilities for qq — > meson and qqq — > baryon coalescence. 
The coefficients C can in general depend on momentum, for example, if there is a 
strong radial flow. 

Equation (JTJ is valid for rare processes only. At high constituent phase space 
densities, most quarks recombine into hadrons and the number of hadrons is linearly 
proportional to that of quarks, dNh(p±_) oc dN q (pj_/n) (n = 2, 3 for mesons, baryons). 

At lower constituent densities, parton coalescence is relatively rare and Eq. (JTJ 
is applicable. Clearly, most partons hadronize via independent fragmentation in this 
case. Still, hadron production can be dominated by coalescence in a broad region 
of hadron transverse momenta. The reason is that in independent fragmentation 
hadrons with a given p± come from rare high-momentum partons with momentum 
p±/z (z < 1), i.e., dN[ ra9 (p±) ~ dN q (p±/z), whereas via coalescence they can come 
from more abundant low- momentum partons with momenta p±/n, i.e., dN^ oal (p±) = 
C h [dN q ( P± /n)] n . 

At very high transverse momenta, the fragmentation process wins because due 
to the quadratic/cubic dependence in Eq. 0) the coalescence yield drops steeper as 
a function of p± than the fragmentation yield. E.g., a power law parton spectrum 
dN q / P± d P± ~ Apl a implies dN% oal /dN[ rag ~ C h A n - l p]^ n ~' l)a -> at high p ± . The 
boundary between the coalescence and fragmentation regions may be at rather large 
transverse momenta pj_ ~ 5 GeV 25 . 



3. Elliptic flow amplification and ordering 

In the coalescence region, meson and baryon elliptic flow are given by that of partons 
via[lj 

v 2 m(p±) ~ 2u 2 , g (^), v 2 . B (p±) ~ 3w 2 , g (^) , (2) 

as follows from Eq. JQ and »2 < 1. For example, if partons have only elliptic 
anisotropy, i.e., dN q /p±dp±d$ — (l/2ir)dN q /p±dp±[l + 2w 2 . g cos(2$)], then 

, s 3tt 2 , g (p±/3) + 3i£ 9 (p ± /3) 2f 2;g ( P± /2) 
1 + 6< 9 (p_l/3) l + 2vl g (p±/2) 

Note, Eq. (0) is also valid for all flow coefficients = (cosfc^), not only u 2 . 

Fig. [2^ illustrates the effect of parton coalescence on baryon and meson elliptic 
flow, as given by Eq. J2J), compared to parton elliptic flow. The latter is shown 
schematically by the solid line. At small transverse momenta, parton i> 2 (p_i_) oc p^_, 
as follows from general analyticity considerations. This region, before v 2 becomes 
approximately linear in p^ could be relatively small (depending on the effective mass 
of partons). At higher transverse momenta p± > 1 — 2 GeV, parton elliptic flow 
saturates as predicted by the MPC parton transport model 3 . 

In the low-pj_ region where w 2 (p±) increases faster than linearly, v% b < v 2,m < 
V2, q - Though here Eq. may not be applicable, it is nevertheless interesting that the 
data do exhibit such a behavior. This ordering follows naturally from hydrodynamics, 
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Figure 2. Qualitative behavior of baryon and meson elliptic flow as a function 
of p± from parton coalescence, a) if all partons have the same V2, b) if < V%. 



where flow decreases with increasing particle mass|5I ITTlll(J| . Similar mass dependence 
could also arise in a coalescence model because heavier hadrons can be formed by 
quarks with larger relative momentum (ignored in the local coalescence approach). 

On the other hand, at highp^i where parton v 2 (j>i_) increases slower than linearly, 
baryon flow becomes larger than meson flow, v 2 b > v 2,m > v 2,q, by as much as 
50%. Parton collective flow saturation above p± > 1 — 2 GeV results in saturating 
meson/baryon flow at p± > 2 — 4 GeV that is amplified two/three-fold compared to 
that of partons. Saturation sets in at 50% higher p± for baryons than for mesons. 
Note that any eventual decrease of parton elliptic flow at very high p± (e.g., because 
of parton energy loss) would happen at two to three times larger p± for hadrons. 

The high-pj^ results above strongly differ from those obtained in Ref. [23] • The 
reason is that, unlike Eq. J2J), in Ref. [231 the coalescence of quarks was considered to 
be independent of their relative momenta and therefore hadron elliptic flow at high 
p± was similar to that of a high-pj_ quark. 

If not all partons show the same elliptic flow, further differentiation occurs because 
in that case 

V2,B=abc(p±) ~ U2,o(Px/3) + «2,&(Rl/3) + «2,c(p±/3) 

V2 : M=ab(p±) ~ V 2 ,a(p±/2) + V 2 ,b(j>±/2) ■ (4) 

For example, strange quarks may have a smaller v 2 (p±) than light quarks, at high pj_ 
because heavy quarks are expected to lose less energy in nuclear medium |27|. while at 
low p± due to the mass dependence of hydrodynamic flow. As illustrated in Fig. yp, 
if i;f < 1)%, elliptic flow decreases with increasing relative strangeness content within 
the baryon and meson bands, i.e., > V2 ~ > v 2 > v 2 and v 2 > v 2 > v^- 

The above conclusions are also supported by detailed numerical calculations using 
the MPC parton transport model|2S]. Fig. [S^ shows the elliptic flow results for 
Au + Au with b = 8 fm at y/s = 200 A GeV at RHIC, computed from the azimuthal 
distribution of hadrons obtained via Eq. l[T|). The initial parton spectra, shown in 
Fig. Et>, were computed from leading-order perturbative QCD (GRV98LO, BKK95, 
K = 2.5, Q 2 = j/|_) with a cutoff po = 2 GeV, below which the spectra were continued 
smoothly down to p± = to yield a p^-integrated parton density dN/drj = 2000 at 
T) = 0. The parton cross sections were cr gg ^ gg = 3 mb = (9 / 'A)o~ qg ^ qg = (9/A) 2 a qq ^ qq . 
The preliminary calculation shows rather small difference between the elliptic flow of 
strange and light quarks. Surprisingly, v 2 is slightly larger than w|, however, this may 
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be because flavor changing channels gg <-> gq and gigi — > (UQa were not included. 
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Figure 3. a) Preliminary elliptic flow ordering results from the MPC parton 
transport model 1551 . b) Initial parton spectra for the MPC simulations. 



4. Solution to the opacity puzzle 

While hadronization via Iparton — > In or independent fragmentation approximately 
preserves elliptic flow at high 2 < p± < 6 GeV , parton coalescence enhances V2 two 
times for mesons and three times for baryons. Hence, the same hadron elliptic flow 
can be reached from 2 — 3 times smaller parton V2, i.e., with smaller parton densities 
and/or cross sections. 

To determine the reduction of the parton opacity quantitatively, one can utilize 
the results of Ref. that computed gluon elliptic flow as a function of the transport 
opacity, X= J dzcr tr p(z) ss adN/drj / (940 mb), from elastic parton transport theory 
for a minijet scenario of Au+Au at y^s = 130A GeV at RHIC. Those results can be 
conveniently parameterized as V2(p±,x) — v ™ ax (x) tanh[p_i_/po(x)]> where u™ aa: is the 
saturation value of elliptic flow, while po is the pj_ scale above which saturation sets 
in. For the estimates here one may assume that all gluons convert, e.g., via gg — ► gq, 
to quarks of similar p± and hence v\ (pj_ ) = v\ (p±_ ) ■ 

Because the increase of elliptic flow with opacity is weaker than linear [31 2], 
v max ^ ^-0.61^ fog charged particle v% data from RHIC can be reproduced in the 
coalescence scenario with 3 — 6 times smaller parton opacities adN/dr](b — 0) ~ 
7000 — 15000 mb than those found in Ref. 3 . The lower(upper) value applies if 
high-pj^ hadrons are mostly baryons (mesons). Preliminary PHENIX data [2H| show 
TTo/h ± ~ 0.5 between 2 < p± < 9 GeV, therefore one may expect mesons /baryons « 
1, in which case adN/dn{b = 0) w 10 000 mb. 

In Ref. |3] only collective flow was considered and the parton opacity at RHIC 
was extracted using elliptic flow data from the reaction plane analysis. Taking into 
account non-flow effects that contributed up to 15-20% [21] to the first elliptic flow 
measurements, parton opacities should be further reduced by 25% to adN/dr](b = 
0) ~ 5000 — 10000 mb. For a typical elastic gg — » gg cross section of 3 mb, this 
corresponds to an initial parton density dN/dn(b — 0) ~ 1500 — 3000, only 1.5 — 3 
times above the EKRT perturbative estimate |2*U]. 
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The remaining much smaller discrepancy is within theoretical uncertainties. For 
example, perturbative cross section and parton density estimates may be too low. If 
most hadrons form via coalescence, the observed hadron multiplicity dN^/dr] ss 1000 
would imply higher initial parton densities dN/dn ~ 2000 — 3000. Constituent quark 
cross sections, a qq w 4 — 5 mb, also point above the « 3 mb perturbative estimate. One 
effect that estimate ignores is the enhancement of parton cross sections a oc a^/fj 2 
due to the decrease of the self-consistent Debye screening mass fi ~ gT e ff(r) during 
the expansion. Finally, the contribution of inelastic processes, such as gg <-> ggg, to 
the opacity has also been neglected so far. A preliminary study shows Q7j that this 
contribution can be similar to that of elastic processes. 
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